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Abstract. The exact 1+3 covariant dynamical fluid equations for a multi-component 
plasma, together with Maxwell’s equations are presented in such a way as to make them 
suitable for a gauge-invariant analysis of linear density and velocity perturbations of 
the Friedmann-Robertson-Walker model. In the case where the matter is described 
by a two component plasma where thermal effects are neglected, a mode representing 
high-frequency plasma oscillations is found in addition to the standard growing and 
decaying gravitational instability picture. Further applications of these equations are 
also discussed. 
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1. Introduction 

Plasmas and electromagnetic fields have an established widespread presence in the 
universe and are known to play an important role in many astrophysical and cosmological 
processes. Although in most cases plasma physics can be adequately addressed within 
the Newtonian or the special relativistic framework, there are occasions where general 
relativistic considerations should be taken into account. The physics of the early universe 
offer a very good example in this respect. General relativistic treatments require the 
rigorous setup of a self-consistent set of equations to describe the plasma dynamics. 
Moreover, when perturbative techniques are employed, there are extra considerations, 
such as those related to the gauge invariance of the approach. In this paper, we will 
try to provide such a setup in the context of cosmological fluid dynamics, leaving the 
possibility of a kinetic-theory based description open for the future.$ 

f When analysing the CMB spectrum, the kinetic approach is used for the photons, while the electrons 
are treated as a fluid (their interaction is mediated via Thomson scattering). This is in contrast to 
many Newtonian applications of plasma physics, where the particle nature of the electromagnetic field 
is neglected, while electrons are described using a kinetic treatment. 
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A number of techniques can be used to analyse the equations describing 
general-relativistic plasmas. Depending on the nature of the problem one might 
employ analytical, numerical and/or perturbative methods. Analytical results are 
usually based on severe symmetry assumptions, which unavoidably restricts their 
applicability. Moreover, the inherent nonlinearity of Einstein’s theory means that 
numerical techniques are also non-trivial to apply. Thus, in many cases the most useful 
method is the perturbative one, possibly combined with numerical methods. In general, 
we may distinguish between two types of approach: 

Non-gravitating plasmas on curved background spacetimes: This method is probably 
best applied to astrophysical situations, and effectively it comprises two sub-cases: (a) 
weak gravitational helds, described by a single potential, or weak gravitational waves; (b) 
strong gravitational helds, where one uses exact solutions to Einstein’s held equations for 
the background. The “membrane paradigm” (see [1]) is a good example of a formalism 
which has been developed for this purpose. 

Self-gravitating plasmas: In this case one takes into account the plasma contribution 
to the total gravitational held. This approach, which is more technically demanding 
than cases (a) and (b) above, is applicable to early universe studies, when most of the 
baryonic matter was ionised. Below, we will give some examples of studies that have 
been based on the above described techniques. 

Considerable amount of work has been done on the interaction between plasmas and 
gravity waves and on the use of electromagnetic helds for the detection of gravitational 
waves (see [2-4] and references therein). This includes studies of parametric excitation 
of plasma waves in the presence of gravitational radiation [5], the scattering of gravity 
waves on highly energetic plasmas during supernovae explosions [6], and the possible 
existence of radio waves due to the emission of weak gravitational waves from binary 
pulsars [7]. Also, in analogy to the frequency upshifting of short laser pulses observed 
in laboratory plasmas (e.g. see [8]), it was shown that weak gravitational waves could 
induce similar phenomena in magnetized multi-component plasmas [9]. Moreover, in [10] 
the exact plane-fronted parallel (pp) solution to Einstein’s held equations (e.g. see [11]) 
was employed to gain a better understanding of nonlinearities in the interaction between 
plasmas and gravitational waves (see also [12]). 

A number of papers employ the membrane paradigm [1], together with the 
appropriate huid equations, to look into the plasma properties in the vicinity of compact 
astrophysical objects such as black holes. In [13], for example, the authors studied 
high frequency EM-waves in plasma outside a spherically symmetric black hole, and 
in [14] they show the possibility of an EM-wave outburst from black holes due to mode 
conversion. Studies looking at the plasma behaviour near rotating black holes can also 
be found in the literature [15]. 

Work has also been done on huid dynamics and kinetic gas theory with the context 
of cosmology. Notably, the book by Bernstein [16], which treats the gas kinetics in 
the Friedmann-Lemaitre-Robertson-Walker (FLRW) model. Nevertheless, there are 
relatively few relativistic cosmological studies that take into account plasma ehects and 
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the behaviour of matter in the presence of electromagnetic helds [17-23]. Thus, the 
general relativistic treatment of plasmas, both in astrophysics as well as in cosmology, 
looks like a field open to investigation. 

When studying relativistic cosmological perturbations, Bardeen’s gauge-invariant 
formalism is the most influential approach [24]. However, Bardeen’s theory is one of some 
complexity and his variables do not always have a transparent physical and geometrical 
interpretation. Moreover, the approach is limited to linear perturbations around a 
FLRW background. Building on the work of Hawking [25] and Olson [26] and utilising 
that of Steward and Walker [27], Ellis and Bruni [28] introduced a mathematically 
simple and physically transparent perturbation scheme. Their formalism, which is both 
covariant and gauge-invariant, has the additional advantage of not been conhned to 
perturbed FLRW universes (see [29] for a comprehensive review). The single fluid 
analysis of Ellis and Bruni has been extended to multi-component systems by Dunsby, 
Brnni and Ellis [31], where a number of possible cosmological applications was discussed. 
Here, we will apply the multi-component formalism of [31] to the case of a charged two- 
component flnid. 


2. Preliminaries 


2.1. The multi-component fluid 


We assnme a family of fundamental observers moving with 4-velocity and a collection 

of perfect flnids with individual 4-velocities given by 

<i) =7(i)(M“ + 'i'p)) , (1) 

where yp) = (1 — is the Lorentz-boost factor and v^^-^Ua = 0 (i is numbering 

each flnid). By assumption each fluid has, in its own rest frame, an energy momentum 
tensor of the form 


^(“) = iTii) + P(i))wp)Wp) + , (2) 

where /i(j) and p(j) are the flnid’s energy density and pressure respectively, while Qab is 
the spacetime metric. Note that in general each species has its own equation of state. 
Relative to the fnndamental frame however, the above reads 


which is the stress-energy tensor of an imperfect flnid with 


( 3 ) 


fl{i) — 7(i)(h'(*) T Pii)) P{i) 1 (4®) 

P{i) = P(i) + llfi) (hp) + P(i))vfi) , (46) 

Cb = 7p)(hp)+Pp))t'p) , (4c) 

= 7p)(hp) + P(b)Ki)4) “ ’ (4^^) 


and u°'u^ is the projection tensor orthogonal to Note that is the heat 

flow and is the anisotropic pressure of each fluid component relative to Clearly, 
both quantities depend entirely on the motion of the species relative to 
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Charged fluids will interact with each other in the presence of an electromagnetic held. 
Thus, we also assume the presence of an electromagnetic held described by the Faraday 
tensor 

, (5) 

where and 5“ = are respectively the electric and magnetic fields 

as measured by the fundamental observers {cabc is the spatial permutation tensor). The 
electromagnetic held contributes to the total energy momentum tensor by 

= \{E^ + B‘^)u^u^ + 1(^2 + , ( 6 ) 

where angled brackets indicate the projected, symmetric and trace-free part of spacelike 
vectors and tensors. Finally, the held obeys Maxwell’s equations 

= (7a) 

V[aF),c]=0. (76) 


2.3. The gravitational field 


The dynamics of the gravitational held is determined by Einstein’s equations, forming 
a closed system once the equation of state for the individual fluid components has been 
established. Of course, in the presence of other physical fields (e.g. anisotropic stresses 
or spinor fields) we need to supplement the system with the corresponding evolution 
and constraint equations (e.g. see [32,33] and references therein). In the presence of an 
electromagnetic held, the conservation laws for the individual charged species are 


= -E%3l^ + 
^0 


1 




( 8 ) 


with = Pc{i)U^i) being the 4-current, pc(i) = the charge density in the rest 

frame of the fluid and Cq is the permittivity of vacuum. The term represents 
interactions other than electromagnetic between the fluids and splits as 


"^ 3 ) = + /(b ) (9) 

where £(*) is the work per unit volume due to the interaction and is the force density 
orthogonal to Because of overall energy-momentum conservation we require that 
J(i) = 0 and write the total fluid equations as 

= (10) 

I I 

Moreover, particle conservation ensures that 

Va(n(qu“i)) = 0 , (11) 

where n(i) is the number density of the individual species in their own rest frame. Finally, 
we point out that the current density in Eq. (8) can be written where 

g(j) is the individual charge of the particles that make up the fluid.§ 

§ In general, we need to employ the second law of thermodynamics '\7aS°‘ > 0, supply an equation of 
state for the species and use the covariant equations given in the Appendix. 
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3.1. The nonlinear equations 

The conservation laws of the individnal flnid components, relative to the frame, are 
obtained by inserting decompositions (4a)-(4o?) into Eq. ( 8 ). In particular, by projecting 
( 8 ) onto we arrive at the energy density conservation equation 

fl(i) = - (/i(i) +P(i)) (^0 + VaUp)) - im+Pii)) (tp) + + n[^pVa7p)) 

-^'p)Va/i(^) + 7(“)^£p) • (12) 

On the other hand, we derive the momentum density conservation equation 
{Ti^) + Pi^)) = -7(“j'V“pp) - 10 (/ip) + pp)) 

— (pp) +Pp)) ('^p)^fe'^p) + 

+7^ (/^p) +pp)) (^'Pl^p) +^'p)^w'^'>7p)) 

~'^p)'^p)^bPp) + 7p)Vcp)(-E'“ + + 7p)Vp)) (13) 

by projecting ( 8 ) orthogonal to Furthermore, the particle number conservation, 
expressed by Eq. (11), takes the form 

hp) = - 0 np) - np)han[*) - 7 (p^ 7 p)np) + Va (7p)«p)'i'p)) • (14) 

Similarly, the total fluid equations (see Eq. (10)) provide the total energy density 
conservation, 

fi = -0(p + p) - Vag“ - 2 hag“ - aV a , (15) 

and the total momentum density conservation 

(/i + p)h“ = -V> - |0g“ - - a%q^ - 

+ (16) 

where pc = X)jPc(i), = 'Thidfi) total charge and current density respectively. 

Also, p = p = Y^iPip, <?“ = Yidp), = Yi^{i) by dehnition and the hatted 

quantities are given by (4a)-(4d). 

The covariant form of Maxwell’s equations is obtained by substituting the Faraday 
tensor, given by (5), into Eqs. (7a) and (7b). They comprise a set of two propagation 
and two constraint equations given by [19-22] 

fl(a} ^ _ | 0 ^a ^ ^ e'^^^oJbE, + + curl5“ - , (17a) 

eo 

5 (a> _ - - curlE“ , (176) 

^aE^ = — Pc + 2u}aB°' , (17c) 

Co 

VaB^=-2UaE\ (17d) 

where curli?“ = e°'^^'VbBc, and analogously for curli7“. 
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We will now linearize the equations of the previous section about a FLRW model. 
Given the homogeneity and the isotropy of the FLRW spacetime, all spatial gradients 
and velocity components orthogonal to must vanish in the background. This implies 
that spatial inhomogeneities are hrst order quantities and that yp) = 1 to hrst order. 
In addition, the symmetries of the FLRW background require that the electromagnetic 
held vanish to zero order as well. This in turn implies, through Eq. (17c), that pc has 
zero background value. Similarly, the shear vorticity u!°' and the acceleration 


also vanish to zeroth order. As a result, Eqs. (12)-(14) linearise as follows 

P(i) = — ^0 + Vanp)j {p{i) + P(i)) , (18a) 

im + P(i}) + '^p)) = -V>(q - (/ip) + pp)) n))) + , (186) 

h(i) = - (^0 + VaCp)) np), (18 c) 

where we have ignored non-electromagnetic interactions between the species (i.e. ep) = 
0 = /))^). Similarly, the total huid equations (15) and (16) reduce to 

/i =-0(/i+ p) - Vag“ , (19a) 

(/i + p)h“ =-V>-|0g“-g“. (196) 

Finally, Maxwell’s equations give 

E“ = -|0E“ + curlR“--j“, (20a) 

Co 

=-|0R“ - curlE“ , (206) 

VaE“ = —Pc , (20c) 

Va5“ = 0. (20 d) 

In many applications, it has proved advantageous to adopt the energy frame, dehned 
by the vanishing of the energy hux,|| 

9“ = ZC.)=')' (21) 

i 

In this frame Eq. (16) reduces to 

(h + p)«“ = -V>, (22) 


which means that for a dust background the acceleration vanishes to hrst order. 

4. Applications 

4 . 1 . Electrically induced velocity perturbations 

Consider an Einstein-de Sitter background and a two-huid system, with each component 
having a dust-like energy-momentum tensor relative to its own frame. In the 

II The electromagnetic contribution to the total heat flux through the Poynting vector e'^^^EbBc vanishes 
to hrst order. 
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background, the only non-zero scalars are the total density /i = and the expansion 

0. Note that to zero order the total charge vanishes (i.e. pc = —e{ni — n 2 ) = 0), since 
both species have equal but opposite charges qi = —e = —q 2 - It follows that pc is a first 
order gauge-invariant variable [27]. Furthermore, pi = mini since no thermal effects are 
included. In this environment, it is useful to introduce the variables 

N = ni+n2, n = ni-n2 , = ^{v'^ + V 2 ), v'" = - V 2 ). (23) 

Given our frame choice (i.e. q°‘ = 0), Eq. (4c) leads to the first order result /iin“ = —^ 2^2 
and subsequently to the following relation 

, ( 21 ) 
p 


between and where 5p = pi — p 2 and 5p/p is independent of time. Then, 
employing Eqs. (186) and (18 c) we obtain the propagation formulae for iV, n and 


N 

h 


- (0 + VaE“) N , 

- 0n - NVaV^ , 

2 mim2 




(22 a) 
(226) 
(22 c) 


As expected, Eqs. (22 a) and (226) show how velocity perturbations, depending on the 
sign of their 3-divergence, can increase or decrease the number density dilution caused 
by the expansion. More importantly, Eq. (22 c) shows that the presence of the electric 
field acts as a source of linear velocity perturbations in the charged plasma, even when 
such perturbations are originally absent (i.e. when Ua = 0 initially). In what follows we 
will see that a non-zero initial velocity perturbation can give rise to density fluctuations 
(cf. (226)), which through Eq. (20 c) may seed electric fields. 


4-2. Velocity induced density perturbations 

Consider the dimensionless, first-order, gauge-invariant variable 

A = ^V^N, (23) 

where a is the background scale factor and flVi, is the covariant Laplacian 

operator normal to u°‘. The above describes inhomogeneities in the total number density 
of the particles and, consequently, it also describes inhomogeneities in the total energy 
density. To linear order the evolution of A is determined by the system 


A = 


(24a) 


p 


= 

- 1^2^-\N 

(mi -I- m 2 )A -I- (mi — m2)a^V^'F 

, (246) 

r = 

- + \a^paY , 

(24c) 

Y = 

--r , 


(24 d) 
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where = 4e^/3eomim2. In deriving the above we have employed the first order 
gauge-invariant variables 

^ = r = aVan“, Y = n/N, (25) 

and used Maxwell’s equation (17c). Note that ^ and Y describe scalar inhomogeneities 
in the expansion and the relative velocity of the species respectively, while Y determines 
the net charge of the total fluid. Given that Eqs. (24c) and (24c?) have decoupled from 
the rest of the system we can obtain the following propagation equation for Y : 

Y + let + la^giY = 0 , (26) 


The solution to Eq. (26) will act as an inhomogeneous driving term in the corresponding 
propagation equation for A: 

A + |0A - l/iA = (|a2 ^ 1) , (27) 

obtained by taking the derivative of Eq. (24a) and using (246). According to Eqs. (24c?) 
and (27), velocity inhomogeneities act as sources of density fluctuations. Note that the 
right hand side of (27) is a pure multifluid effect, where the part containing stems 
from the plasma description. 

In order to solve equations (26) and (27) it is standard to decompose the physical 
(perturbed) fields into a spatial and temporal part, using as eigenfunctions Q{k)i 
solutions of the scalar Helmholtz equation [34]. In particular we write 


A = AnjO'**, r = , 


(28) 


where taY(k) = 0 = Va^{k), Q{k) = 0 and = —{lY /Eor an Einstein-de 

Sitter background, the expansion and energy density evolve as 0 = 2/t and p = A/3C. 
Hence, applying the harmonic splitting given above, Eqs. (27) and (26) become 

~ ~ 

and 

4 • cfl 

tk) + + flitk) = 0 > (30) 

respectively. In order to estimate the value of the parameter a we substitute back for 
the gravitational constant and write 


^2 _ 4 / 1 \ ^ f 

3 \mi) \m2) Veo/ \87rGml) \mi) 
Since a 3> 1 the solutions to above equations are 

A(fc)=Cir2/=' + C2r-i + A;2^F(fc), 

p 

_ 1 

Y(k) = [Cl cos(q; In r)-|-C 2 sin(Q; Inr)] r 6 , 



X 10^^ . 


(31) 


(32 a) 
(326) 


where we have introduced the dimensionless time-coordinate r = t/tj, with ti 
corresponding to some arbitrary initial time. Hence, in addition to the usual growing 
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and decaying modes of the standard gravitational instability picture, we have obtained 
a mode representing high frequency plasma oscillations with a weak damping envelope. 
This mode is triggered by velocity distortions in the charged plasma and, as expected, 
has negligible large scale effect. However, the extra plasma modes become increasingly 
important as we move on to progressively smaller scales (i.e. for /c 3> 1 ). 

It should also be pointed out that a finite temperature will in general cause Landau 
damping of the plasma oscillations. The effect (requiring kinetic treatment) is small for 
wavelengths much larger than the Debye length (which is proportional to the thermal 
velocity of the plasma particles) and in this case the dust fluid approximation is well 
justihed. 

4 . 3 . Velocity induced electromagnetic fields 

In this section, we will derive the wave equations for the electromagnetic field, seeded 
by velocity perturbations. 

For a cold plasma, the currents for each fluid species may be written as 

Jh = + <i)) > (33) 

where g(j) is the charge and v'^-^ is the velocity of the species under consideration. Since 
we require the plasma to be neutral on the whole, the species are of opposite charge. 
Hence, the total current appearing in Maxwell’s equations reads to first order 

f = jt + J2 =-eNv^. (34) 

^From Maxwell’s equations (20a)-(20d), using (34) and (22c), one can then deduce wave 
equations for the induced electromagnetic helds: 

+ |0^(a) + [f©' + (fcT^ + i) /x] Ea = (vaY - lOVa) , (35a) 

B(^a) - + |05(a) + [f©^ + Ba = -2/3^ curl, (356) 

where = e/eo{mi + m 2 ). These equations govern the interaction of density/velocity 
perturbations with electromagnetic waves in the plasma and show in particular that 
density/velocity perturbations induce wave phenomena. Observe that Ba and curl Va 
are both purely solenoidal, whereas VaH has no solenoidal part. It is worthwhile 
to note that the magnetic held is solely sourced by inhomogeneities in the velocity 
in contrast to the electric held which is sourced by inhomogeneities in the number 
density and velocity perturbations. Clearly, the velocity perturbation is non-zero even 
if = 0, as long as 7 ^ 0 initially (cf. (22c)). Both Eqs. (35a) and (356) look strikingly 
similar, the diherences originating either from the total current or from a gradient in 
the charge density (in the case of VaH). The additional 3Q;^/4-term in the electric wave 
equation comes from the non-stationarity of the total current, and its large magnitude 
— ~ 10 ^^ for an e’''e“-plasma — leads directly to the high-frequency behaviour of 

plasma ehects, as will be shown below (see also the preceding section). 
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It will be useful to introduce expansion normalized variables, 


0 ’ 


_ Ba ^ _ curlua 
“ 0 ’ 0 


(36) 


Equations (35 a) and (356), together with equations for the driving terms, then 
4> - VV, + (0 - t) 4) + [-|0' + (|«2 + 1) p] = 2(3%{VaY - \Qva) , 
h(„) + iQVa = -|f^0<^a , 

4a> - +(©-#) 4a> + hi©' + 4] = -2/9^ , 

+ (jS - if) = If + (|e - 5«)« " 


read 

(37a) 

(37b) 

(37c) 

(37(i) 


Equation (37d) follows from (36) using (376) and Maxwell’s equation (206). 

Restricting ourselves to scalar perturbations, we take the divergence of the above 
equations to extract the scalar part of the system. Of course, there is no contribution 
from the magnetic held in this case. Using (25) and dehning S’ = aV“(^a, Eq. (376) then 
transforms into (cf. (24c)) 


r+i0r =-|^0^= , (38) 

where the last equality is a direct consequence of Maxwell’s equation (20c). Combining 
Eq. (24d) with (38) and using (19 a) together with the commutator expression 


aV“VVa = VV + (-§02 + 4) ^ , 


(39) 


one can show that the scalar part of the electric wave Eq. (37a) reduces to^ 

i + (40 _ I ) ^ + [§02 + (|«2 _ 1 ) ^ ^ 0 _ ( 40 ) 

In addition, Eq. (38) gives rise to propagation equations for Y and Y, as discussed 
earlier: 


r+lQy+[-lQ^+{la^-l)fi\y = 0, (41a) 

U +f0y + = 0 . (416) 

Hence, Eqs. (40)-(416) all stem from (38). 

Specialising to a hat FLRW model with a zero cosmological constant, for which 
/i = 1/3 0^ and 0 = 2/t always holds, solutions to these equations can easily be 
obtained: 


Y{t) = ^ {Acos{ijj\nT) + ^ {^\A + B) sm.{ijj\nT)^ , (42a) 

= -\^^{{‘^^ + ^B)cos{(3:\nT) + sin(o;lnr)| , (426) 

{(2^ + 3H)cos(o;lnr) + sin(o;lnr)| . (42c) 

Here, we used again the dimensionless time-coordinate r = t/U, where ti denotes some 
arbitrary initial time. Initial conditions of the velocity perturbation are chosen to be 


^ Note that in deriving Eq. (40), the Laplacian terms cancel, and a harmonic decomposition is therefore 
not needed. Thus, the electric field will not contain a particular length scale, due to its Coulomb-like 
nature. 
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A = and B = Yfll) (a prime stands for dfl)- The frequency of the solutions 

is proportional to uj = and grows logarithmically in time. The solutions 

exhibit high-frequency plasma behaviour. Observe that although the solutions decay 
with time, their magnitude changes only very slowly over time, particularly if the velocity 
perturbations are taken to start at the onset of recombination. 

We have restricted our attention to scalar perturbations, with the implication that 
magnetic held effects vanish. From the point of view of generating magnetic seed helds, 
for, e.g., the dynamo mechanism or the Biermann-battery effect (see [30] and references 
therein), it is of interest to analyse vector perturbations in a similar way. This is reserved 
for future research. 

5. Discussion 

In this paper we generalized the multi-component huid equations derived by Dunsby, 
Bruni & Ellis [31] to the case of charged huids in the presence of electromagnetic helds. 
The equations are given in covariant form, relative to an observer moving with velocity 
that is taken to coincide with the average velocity of the cosmic medium. We linearized 
these equations about a FRW universe and then applied them to an Einstein-de Sitter 
(EdS) universe. Our matter held is an ion-electron plasma with zero average pressure 
(which made the EdS model a suitable background). We showed how, when there is a 
residual net charge, the presence of an electric held can lead to velocity perturbations 
even when the latter are originally absent. We also found that velocity distortions can 
source inhomogeneities in the number density, and therefore in the energy density, of the 
huid. In fact, our linear equations reveal the presence of an extra mode, representing 
high frequency plasma oscillations, in addition to the standard growing and decaying 
modes. This mode is likely to be important on scales considerably smaller than the 
Hubble radius and therefore is of little importance as far as structure formation is 
concerned. It does illustrate, however, interesting small scale physics that could play a 
role during the latter stages of galaxy formation. 

We also applied our covariant equations to look into the generation of 
electromagnetic helds due to velocity perturbations in a plasma. The corresponding 
wave equations, with the velocity distortions playing the role of a source, were given, and 
they were solved in the case of scalar perturbations. The solutions show high-frequency 
behaviour typical of a plasma. We restricted our attention to scalar perturbations, thus 
obtaining an evolution equation for the electric held. However, magnetic held ehects 
were absent since these are related to vector modes. Because magnetic seed helds play 
a crucial role in, e.g., the dynamo mechanism, it is of great interest to pursue the 
analysis of the presented equations in the context of vector perturbations. Results in 
this direction will be presented elsewhere. 

There are a number of ways to generalize the discussion presented in this paper. One 
possibility is to include thermal ehects which occur in a photon-baryon plasma giving a 
non-zero acceleration to hrst order. This may lead to possible coupling between acoustic 
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and plasma oscillations. In addition one could apply the ponderomotive force concept 
between neutrinos and electrons (see [32,35] and references therein) to cosmology in a 
covariant context. In this picture, derived from the theory of electroweak interactions, 
there is an effective interaction between electrons and neutrinos due to density gradients 
in either species. For instance, the (non-relativistic) force density exerted by neutrinos 
on the electrons is given by [35] 

fC = —^ (l + 4sin^ Ow) GprieV^'n^ , (43) 

V2 

where 6 w is the Weinberg angle and Gp is the Fermi constant. The expression 
(43), together with its neutrino counterpart, could act as a driving force for density 
fluctuations in the early Universe, possibly giving a neutrino signature in the CMB, 
having an alternating structure as compared to the regular CMB spectrum. The 
neutrino-driven instability discussed by Silva et ah [36] (see also Ref. [37] for the 
covariant relativistic form of the same equations), using kinetic theory, could in principle 
be transferred to a gauge invariant covariant formalism, suitable for cosmological 
applications (see also [38]), but this is left for future studies. 
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Appendix A. Gravitational dynamics 

The covariant equations for the dynamics of the gravitational held was given in Ref. [29], 
and we use their notation. 


Appendix A. 1. Covariant equation 

• Evolution equations for kinematic variables: 

0 — V^ 0^ + — 2 -|- 2 ufl — ^ (/^ + 3p) -|- A , (-^-1) 

^(a> _ I ^abc | 0 ^ (A.2) 

(jW _ fljAyb) = _ I 0 + uC yb) _ ^b)c _ ^{a ^b) _ ^^ab _ 1 ^ab^ ^ ^^.3) 

• Constraint equations for kinematic variables: 

0 = - I V“0 + [ Vboo, + 2u^u,]+q\ (A.4) 

0 = - (h,o^“) , (A.5) 

0 = 77 “^ + 2 - (curl , (A.6) 

where (curl (t)“^ = 
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• Evolution equations for the curvature variables: 

- ( curl Hfl’’ + i = - 1 (^ + p) - 0 + i 

+ 3 i [ 2 ii, H% + o;, + \ A) ] , (A.7) 

+ (curl - i( curl - 0 + 3 + | 

- [ 2 ii, E% - I q, - u, H% ] , (A.8) 

where 

{curl Hfl’’ = 7]'^'^^^ VcH’^^d , (A.9) 

( curl E)“^ = , (A. 10) 

( curl 71 VcA • (A. 11) 

• Constraint equations for the curvature variables: 

0 = Vb(E“'' + 17r“'’) - i V> + i 0 - I g'’ - 3 

- [ ^bd - I cufe gc ] , (A. 12) 

0 = + (/i + p) o;“ + 3 cufe - 17r“') 

+ 7]°'^'^ [ I Vfegc + Ufed (E'^c + I '?i''c) ] • (A. 13) 
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